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Abstract. We present new spectroscopic and photomet- 
ric time series observations of the 5 Scuti star FG Vir. We 
detect the oscillations via changes in the equivalent widths 
of hydrogen and metal absorption lines. From the ratios 
between spectroscopic and photometric amplitudes, we as- 
sign £ values to the eight strongest oscillation modes. In 
particular, we identify two radial modes (£ = 0) and find 
that the main pulsation mode (147 /d3z) has £ = 1. One of 
the radial modes (at 140/iHz) is the fundamental, imply- 
ing that two modes with lower frequencies are g-modes. 
For the radial modes, we compare frequencies with those 
calculated from a scaled S Scuti star model and derive a 
density 0.1645 ± 0.0005 p©. We then obtain a distance of 
84±3pc, in excellent agreement with the Hipparcos value. 
Finally, we suggest that a 3.5-day variability in all observ- 
ables (equivalent widths and intensity) is caused by stellar 
rotation. 

Key words: stars: S Scuti - stars: oscillations - stars: 
individual:HD106384=FG Vir - techniques: spectroscopic 



1. Introduction 

Oscillations in multi-periodic variables such as 6 Scuti, 
roAp and (3 Cephei stars have been observed extensively 
during the past 20 years. But even with high-quality data, 
it is still extremely difficult to identify which modes are 
being detected. Kjeldsen et al. ( 1995| ) used a new tech- 
nique to detect solar-like oscillations in the bright G sub- 
giant r\ Boo through their effect on the equivalent widths 
of the B aim er lines. A subsequent discussion by Bedding 
et al. ( 1996 ) of the sensitivity of different observables to 
modes with different degree £ suggested that one can de- 
termine the ^-value of a given mode by combining mea- 
surements of absorption-line equivalent widths with simul- 
taneous photometric observations. To test this idea, we 
chose the bright and well-studied 5 Scuti star FG Vir. 

FG Vir (HD 106384; V = 6.57) is a multi-periodic 
6 Scuti star. It has a main pulsation period close to 1.9 



hours and shows a fairly complex oscillation spectrum. 
This star has been studied extensively during the last 
few years, resulting in the detection of at least 24 well- 
determined frequencies between 100 and 400 pTTz, with 
amplit udes fr om .8 to 22 milli-magnitudes (mmag; Breger 
et al. 



1995 



Because of its slow rotation (which 
reduces the complicating effects of rotational splitting) 
and the large number of detected frequencies (some of 
which are probably g modes), FG Vir is one of the most 
promising candidates for performing asteroseismology on 
a 5 Scuti variable. Observations and m odels of this star 
hav e bee n presented by Dawson et al . ( 1995 ), Breger et 
al. (|1995D and Guzik & Bradley fll995D. 



By choosing a star like FG Vir we have the advantage 
of knowing the frequencies in advance. We are therefore 
able to determine the oscillation amplitudes and phases 
with high precision. The aim of this paper is to identify 
the £ values of the observed modes and to compare the 
oscillation frequencies with a pulsation model. A prelimi- 
nary analysis of the observations presented in this paper 



was given by Viskum et al. (1997b), while results on ra- 
dial velocity measurements were given by Viskum et al. 
(|l99"7al). 
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2. Observations and basic data reduction 

2.1. Spectroscopy 

We obtained intermediate-resolution spectra of FG Vir 
using the DFOSC spectrograph mounted on the Danish 
1.54-m at La Silla, Chile and the coude spectrograph (B 
grating) on the 74-inch Telescope at Mt. Stromlo, Aus- 
tralia, during a period of more than two weeks in May 
1996. A journal of the observations is given in Table 0. In 
total we obtained 652 spectra of FG Vir. 

The DFOSC spectrograph was used in the echcllc 
mode, covering a wavelength region from 4600A to 8000A 
at a resolution of R ~ 4300. Six echelle orders were 
recorded, each covering about 600A and with an over- 
lap between adjacent orders of about 370 A. A S/N ratio 
of 250-360 was reached in the DFOSC data with a typi- 
cal exposure time of 460 s. The S/N ratio was determined 
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Table 1. Journal of spectroscopic observations of FG Vir Table 2. Journal of photometric observations of FG Vir 



Date 




Site 


HJD (start) 


Length 


No. of 


Date 




HJD (start) 


Filter 


Length 


No. of 


UT 






2450000 + 


hours 


spectra 


UT 




2450000 + 




hours 


data pts 


1996 May 


2 


Stromlo 


205.931 


4.0 


26 


1996 Apr 


28 


202.469 


uvby-/3 


7.1 


62 


1996 May 


3 


Stromlo 


206.834 


4.2 


31 


1996 May 


1 


205.470 


uvby-/3 


6.9 


69 


1996 May 


6 


Stromlo 


209.867 


7.7 


62 


1996 May 


2 


206.472 


uvby-/3 


6.8 


46 


1996 May 


7 


Stromlo 


210.838 


5.9 


47 


1996 May 


3 


207.465 


uvby-/3 


7.2 


48 


1996 May 


8 


La Silla 


212.484 


6.2 


34 


1996 May 


4 


208.490 


uvby-/3 


6.4 


42 


1996 May 


9 


La Silla 


213.486 


6.3 


45 


1996 May 


5 


209.469 


Hj3 


7.0 


83 


1996 May 


10 


La Silla 


214.496 


6.1 


43 


1996 May 


7 


210.511 


HP 


5.2 


51 


1996 May 


11 


La Silla 


215.496 


6.0 


50 


1996 May 


8 


212.464 


uvby-/3 


6.9 


46 


1996 May 


13 


La Silla 


216.550 


4.7 


27 


1996 May 


9 


213.465 


uvby-/3 


3.2 


22 


1996 May 


13 


La Silla 


217.492 


6.0 


41 


1996 May 


10 


214.457 


uvby-/9 


7.3 


48 


1996 May 


15 


La Silla 


218.525 


5.1 


40 


1996 May 


11 


215.458 


uvby-/3 


7.4 


49 


1996 May 


15 


Stromlo 


218.824 


6.2 


69 


1996 May 


13 


217.462 


uvby-/3 


7.2 


45 


1996 May 


15 


La Silla 


219.468 


5.1 


35 


1996 May 


14 


218.486 


uvby-/3 


6.4 


42 


1996 May 


17 


La Silla 


220.541 


3.3 


20 


1996 May 


15 


219.470 


H(5 


4.2 


51 


1996 May 


17 


La Silla 


221.477 


6.2 


43 


1996 May 


17 


221.461 


uvby-P 


6.4 


43 


1996 May 


18 


La Silla 


222.498 


5.5 


39 


1996 May 


18 


222.475 


uvby-/3 


3.1 


21 



by measuring the standard deviation of the intensity of 
the normalized spectrum in two spectral windows free of 
absorption lines. 

The Mt. Stromlo data consisted of single-order spectra, 
with a wavelength coverage of 6000A to 7000A. The dis- 
persion was 0.49A/pixel, with a resolution of about 1.5A 
which was set by the slit width of 2 arcsec. The typical 
exposure time was 400 s. 

We used the IRAF package for bias subtraction, flat- 
field correction, extraction of the ID spectra and sub- 
traction of background scattered light. We found the Mt. 
Stromlo CCD to have a non-linearity of 7% from low-light 
levels to saturation, which we corrected in pre-processing. 
The DFOSC CCD camera was measured to be linear 
within 0.5% and no correction for non-linearities was per- 
formed. 

The DFOSC spectra suffered from sky contamination 
and edge vignetting, which meant that the extracted spec- 
tra tapered from the centre to the edge. We therefore nor- 
malized the continuum by dividing each order by a curve 
containing the shape of the flat field. The Mt. Stromlo 
spectra did not suffer for this kind of edge vignetting, so 
a third order polynomial fit was used for the continuum 
normalization. 



2.2. Photometry 

We obtained Stromgren uvby-[3 photoelectric photometry 
with the six-channel photometer mounted at the Danish 
0.5-m reflector (SAT) at La Silla, Chile over 22 nights 
in 1996. Observations of FG Vir were interspersed with 
observations of the star HD 105912 (spectral type F5 V). 
This was used as a comparison star by Breger et al. ( |1995| ), 
who found no variability. In Table [2| we present a journal 
of the observations. The precision per data point is about 



3 mmag. During three of the nights, only the H/3 index 
was obtained due to non-photometric weather conditions. 
A total of 768 useful data points (99 hours) were obtained. 

To produce light curves from the raw data, we used a 
reduction program developed for the SAT telescope. The 
reduction included airmass correction using standard ex- 
tinction coefficients and transformation to the standard 
instrumental system. Slowly varying effects were then re- 
moved by subtracting second-order polynomials for both 
FG Vir and the comparison star. Variations in FG Vir not 
due to oscillations were removed by subtracting the com- 
parison star, resulting in the light curves whose analysis 
is described in Sec. 0. 



3. Determination of equivalent widths 

Stellar acoustic oscillations (p modes) can be measured 
via their effect on the equivalent widths (EWs) of the 
Balmer lines, which are sensitive to temperature. Metal 
lines a re als o sensitive to temperature changes (Bedding 
et al., 1996 ), with the Fel lines being particularly sen- 
sitive but in the opposite sense to the Balmer lines. We 
measured the equivalent widths of the metal lines by fit- 
ting profiles, but this technique was found unsuitable for 
the much stronger Balmer lines, for which we used a dif- 
ferent method. Both techniques are described in the next 
sections. 



3.1. Equivalent-width variations of the Balmer lines 



As described by Bedding & Kjeldsen (1998), a robust way 



of estimating the EW of broad spectral lines involves a 
method analogous to Stromgren H/3 photometry. We have 
attempted several other methods, such as fitting directly 
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to the line profile, without reaching the same internal pre- 
cision. 

Our adopted method involved calculating the flux 
(number of counts) in three software filters, one centred 
on the line (L) and the others on the continuum on the red 
(Cr) and the blue (Cb) sides of the line. For each spec- 
trum, we first adjusted its slope so that the fluxes were 
equal in the continuum filters: Cr = Cb = C. We then 
calculated the quantity 



W = 



(C-L) 
C ' 



(1) 



This calculation was repeated for a large number of po- 
sitions around the line centre by moving the three filters 
together in wavelength by small amounts to find the po- 
sition which maximized the value of W . The whole pro- 
cedure was repeated for different filter widths, in a way 
analogous to aperture radii selection in CCD photome- 
try, and the results combined to give an estimate of the 



equivalent width (see Bedding & Kjeldsen 1998 for more 
details). 

Our La Silla spectra cover both the Ha and H/3 lines, 
while only Ha is present in the Mt. Stromlo data. In 
Fig. |l|a we present the amplitude spectrum of the Ha EW 
for the combined data set (Mt. Stromlo and La Silla), 
where the amplitude is measured in promille (parts per 
thousand). All amplitude spectra in this paper were cal- 
culated using a weighted least-squares sine- wave fit to the 
time series (see Frandsen et al. 1995, 1996| ). The pulsation 
signal is clearly seen in the range 100-400mHz. Fig. [l|b 
shows the residual spectrum obtained after subtracting 23 
previously published frequencies (as described in Sec. ||). 
The mean noise level in the residual spectrum at high fre- 
quencies (600-800/iHz) is 0.18 promille, while it is 0.35 
promille in the frequency range where the oscillation sig- 
nal is found (100 - 400 /zHz). The S/N of the dominant 
peak at 147 /zHz is about 58. 

Fig. H shows the amplitude spectrum of the EW mea- 
surements of the H/3 line (La Silla only). The oscillations 
are again obvious, but the noise is more than a factor two 
higher than for Ha, due in part to bad CCD columns that 
could not be removed by bias subtraction or flat fielding. 



3.2. Equivalent- width variations of metal lines 

The spectra from La Silla consist of six echelle orders con- 
taining many metal lines. The line identifications, kindly 
provided by F. Kupka and M. Gelbmann, were made us- 



ing the program described by Gelbmann (1995) and Gelb- 



mann et al. (1997). We measured the EWs by fitting Gaus- 
sians to the line profiles, with the surrounding continua 
being fitted using low-order polynomials. The fitting was 
done iteratively using the Levenberg-Marquard method 



(Press et al. 1989) to estimate the line position, line depth 
and continuum level. A first estimate was made for these 
parameters over a narrow wavelength interval, and an im- 
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Fig. 1. (a) Amplitude spectrum of the oscillations in 
FG Vir measured through the EW of the Ha line of the 
total data set (Mt. Stromlo and La Silla). The amplitude 
is measured in promille (parts per thousand). The inset 
shows the amplitude spectrum of the window function, 
(b) Residual spectrum after subtraction of 23 known os- 
cillation frequencies. 




200 400 600 

Frequency (/U.Hz) 



800 



Fig. 2. Amplitude spectrum of the oscillations in FG Vir 
measured through the EW of the H/3 line (La Silla data 
only). 
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proved continuum was calculated over a broader range of 
the spectrum. In addition to the EW, this also yielded 
other parameters to be used for decorrelation. 

The time series of each line was decorrelated against 
external parameters. A detailed discussion of our decorre- 
lation procedure can be found in Frandsen et al. ( 199(f ). 
We used parameters that were derived either from the 
spectra or from the CCD frames and were known not to 
contain any information on the EWs or the pulsations. 
Any correlation between the measured EWs and these pa- 
rameters must then have been due to a systematic effect 
and was corrected by subtracting a simple linear fit. The 
procedure was repeated until no significant correlations 
were present. This procedure increased the S/N signifi- 
cantly: for example, the improvement was a factor of 3.3 
for the sum of all 20 lines described below. 

From the EW measurements we detected the strongest 
pulsation period of FG Vir in all Fei, Cai and Mgi lines, 
with amplitudes from 7 and 22 promille. Noise levels in 
the amplitude spectra ranged from about 0.9 promille up 
to more than 3 promille for the weakest or highly blended 
lines. 

We also examined lines of ionized elements, which are 
potentially useful for calibration. The reason is that in FG 
Vir, these lines are near their maximum strength as a func- 
tion of temperature (Gray 1992) and thus are expected to 
be stable with respect to variations in the temperature 
induced by the stellar pulsation. We indeed found these 
lines to have very low amplitudes in EW, with some of 
the Sin lines having no detectable variability. For exam- 
ple, the Sin line at 6372A was stable in EW at the level 
of 2.7 promille (3<r). We used the EW of this line as a 
decorrelation parameter for the Fe I lines. 

The Fei lines give the best opportunity for mode iden- 
tification. Table [| lists the 20 strongest unblended lines, 
where Column 1 is the central wavelength, Column 2 is 
our line ID number (the first digit indicates echelle order 
in which the line falls), Column 3 is the relative line depth 
(in percent), Column 4 is the amplitude of the strongest 
mode and Column 5 is the uncertainty in this amplitude. 
Three of the chosen Fe I lines are detected in two echelle 
orders. For these lines, the two measurements of oscilla- 
tion amplitude agree within the errors. All the lines have 
the same phases within the uncertainties. As evident from 
Table ||, the amplitude of the main mode changes between 
the different lines. This could be due to variations in the 
temperature sensitivity of the different Fei lines, or due 
to blending with other species with different temperature 
sensitivities to Fe I. 

Fig. ^| shows the amplitude spectrum obtained by com- 
bining the EW time series of the twenty strong Fei lines, 
weighted according to their signal-to- noise ratios. To be 



Table 3. Amplitude of the strongest oscillation mode in 
FG Vir as measured in EW for twenty strong Fei lines. 



A 


Line 


depth 


A147 


a(A 147 ) 


(A) 


ID 


(%) 


(promille) 


4958 


1-14 


23 


13.4 


1.1 




2-26 


24 


13.0 


1.5 


4983 


2-25 


14 


18.5 


1.2 


5003 


2-23 


15 


15.1 


1.7 


5041 


1-19 


15 


13.1 


1.2 




2-21 


16 


16.1 


2.1 


5227 


2-13 


25 


12.7 


0.8 


5270 


2-28 


21 


12.4 


0.9 


5329 


2-10 


20 


14.1 


0.9 


5341 


3-25 


20 


13.4 


0.9 


5370 


2-08 


13 


16.3 


1.6 


5404 


3-22 


10 


20.2 


1.8 


5406 


3-21 


11 


21.2 


2.0 


5424 


3-19 


11 


9.6 


2.0 


5430 


3-18 


12 


17.7 


2.3 


5456 


2-04 


11 


6.8 


3.0 




3-15 


13 


9.0 


2.0 


5477* 


2-03 


11 


17.8 


2.2 


5616 


3-05 


11 


15.1 


1.7 


5659 


3-03 


10 


22.0 


1.6 



'this line is blended with a strong Nil line. 
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Fig. 3. Amplitude spectrum of twenty combined Fei lines 
measured through the EW using profile fitting. 

precise, the combined time series w(t) was obtained from 
the 20 individual series Wi (t) using the following equation 



w{t) = 



E^(s/N)?) (SA(S/N); 
(E(s/nx 2 ) 2 



(2) 



where Aj is the amplitude of the 147 /iHz mode and 
S/Nj = Ai/criAi) is its signal-to-noise ratio. The noise 
level in Fig. |3| is 0.5 promille in amplitude. Oscillations 
are present with sufficient signal-to-noise in this combined 
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Table 4. Pulsation modes detected in different spectral 
lines measured via the EW 



ID 


Frequency 




Amplitudes 










Ha line 


H/3 line 


Fei lines 




/iHz 


c/d 


promille 


promille 


promille 


h 


147.18 


12.716 


20.3 


20.2 


14.1 


h 


280.42 


24.228 


3.9 


3.0 


3.1 


h 


270.87 


23.403 


3.9 


3.4 


4.0 


h 


243.65 


21.052 


3.2 


3.3 


1.5 


h 


229.95 


19.868 


3.1 


2.3 


1.4 


h 


140.67 


12.154 


3.7 


3.6 


4.0 


h 


111.76 


9.656 


4.1 


4.3 


1.9 


h 


106.47 


9.199 


3.0 


4.0 


1.6 



data set to allow the amplitudes to be used for mode iden- 
tification. In the next section we measure amplitudes for 
as many modes as possible in each of the observables (EWs 
and photometry) , to allow ratios to be calculated for mode 
identification. 

4. Frequency Analysis 

Our data set suffers from a poor window function, which 
produces very high side- lobes (see Fig. [I]), which compli- 
cates the determination of frequencies in the Fourier do- 
main. Fortunately, Breger et al. ( |1995| , |1998| ) have deter- 
mined very precise oscillation frequencies for FG Vir from 
data obtained over many years. We have therefore chosen 
to rely on these frequencies instead of trying to determine 
them from our data set. 

First, however, we must check whether our data set is 
consistent wi th the main oscillation frequencies found by 
Breger et al. ( 1995 , 1998 ). By pre- whitening with their fre- 
quencies, we foun d tha t the ten strongest pulsation modes 
of Breger et al. ( 1995 ) are also present in our data set. 
From the residual spectrum it was clear that more modes 
are present at lower S/N. 

Knowing the frequencies in advance reduces the task 
of determining amplitudes and phases to a linear problem. 
We used the program ISF (Frandsen et al. 1995 ), which 
is a least-squares routine that simultaneously fits a series 
of known frequencies to a given time series. To maximize 
the precision, we have fitted all 23 frequencies determined 



by Breger et al. (1998) to each time series, except one 
frequency (/20) which is very close to an alias of one of 
the other frequencies. Table [| lists the amplitudes for Ha, 
H/3 and the combined Fe 1 spectral lines. We have used the 



same mode numbering as Breger et al. (1998). Although 
all 23 frequencies were fitted, we shall only use the eight 
modes for which S /N> 5 in both Ha EW and photometry 
(see also Table [5]) . Fig. |l|b shows the residual spectrum for 
Ha EW after subtraction of the fit. 

Fig. ^ shows amplitude spectra (in mmag) for the four 
Stromgren filters. The signal can clearly be seen in each 
filter, with the expected decrease in amplitude from blue 



to red. Eight frequencies are detected with S/N> 5, with 
amplitudes given in Table |. The same eight frequencies 
were detected in the b — y index with lower S/N. In the (3 
index, only two frequencies could be found with adequate 
S/N but the table includes the amplitudes for all eight 
frequencies for completeness. The mean noise levels in the 
residual amplitude spectra of the u, v, b, y, b — y and /3 
time series are 0.42, 0.42, 0.40, 0.35, 0.16 and 0.29 mmag, 
respectively. These were used to calculate the S/N values 
in Table |. 

Because the frequencies are known, the amplitudes can 
be estimated with better precision than would be expected 
from the above noise values. The formal solution gives 
standard deviations of the derived amplitudes for the uvby 
data to be 0.07 mmag. However, from tests made by de- 
termining amplitudes from subsets of the different time 
series, we find that the main uncertainty on the ampli- 
tudes comes from the very complicated window function. 
Furthermore, some of the modes are closely spaced in fre- 
quency and are not well separated in our data. This is es- 
pecially true for the closest pair of modes fi = 280.42/iHz 
and /n = 280.09/iHz, which are unresolved by our ob- 
servations. The derived amplitudes of these modes change 
when data are added to or removed from the time series. 
In order to minimize the influence of the window function 
on the amplitudes, so as to be able to compare amplitudes 
from different observables, we have removed a few nights 
from the photometric time series to obtain an observing 
window similar to that of the La Silla spectroscopic obser- 
vations. In this way, even if amplitudes of some modes are 
slightly wrong due to the complicated window function, 
this effect should be about the same in each time series 
and therefore cancel when we calculate the amplitude ra- 
tios. 

Breger et al. ( 1995] ) presented broad-band V photom- 
etry for the ten strongest modes in FG Vir, which we 
can compare with our Stromgren y amplitudes. For the 
strongest mode (/1) they obtained an amplitude of 22.0 
mmag, which is in good agreement with our y amplitude. 
We also agree on all other amplitudes, within the observ- 
ing errors, except the J/3 mode at 270/iHz, which has a 
much higher amplitude in our data set. As discussed by 



Breger et al. (1998), the amplitude of the f% mode has 
been highly variable. It has increased from 1.4 mmag in 
1992, 2.3 mmag in 1993 to 4.1 mmag in 1995, and our 1996 
measurement gives 5.2 mmag. 

5. Mode identification 



The investigation by Bedding et al. (1996) into the sen- 
sitivity of various observing methods to modes with dif- 
ferent I showed that the Balmer-line EWs have a spa- 
tial response similar to that of radial velocity measure- 
ments. This is because velocity measurements only detect 
the component along the line of sight to the star and so, 
when integrated over the stellar disk, are preferentially 
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Table 5. Stromgren photometry for FG Vir 



ID Freq. Amplitudes 



u v b y b — y (3 





/iHz 


mmag 


S/N 


mmag 


S/N 


mmag 


S/N 


mmag 


S/N 


mmag 


S/N 


mmag 


S/N 


h 


147.18 


32.5 


77 


32.4 


77 


27.6 


69 


22.0 


63 


5.4 


33 


5.3 


18 


h 


280.42 


6.5 


15 


6.5 


15 


6.5 


16 


5.7 


16 


1.2 


7.3 


0.7 


2.5 


h 


270.87 


7.1 


17 


7.3 


17 


6.5 


16 


5.2 


15 


1.0 


6.5 


1.3 


4.4 


h 


243.65 


4.6 


11 


4.5 


11 


4.2 


10 


3.6 


10 


0.7 


4.2 


1.2 


4.0 


h 


229.95 


3.3 


7.8 


3.3 


7.8 


2.8 


7.0 


2.6 


7.5 


0.4 


2.5 


0.5 


1.8 


h 


140.67 


7.8 


18 


7.6 


18 


6.0 


15 


5.3 


15 


0.8 


5.1 


1.0 


3.4 


h 


111.76 


5.3 


12 


5.2 


12 


5.0 


13 


4.3 


12 


0.6 


4.0 


0.9 


3.1 


h 


106.47 


3.9 


9.3 


4.3 


10 


3.5 


8.8 


3.0 


8.6 


0.5 


3.3 


0.8 


2.8 



weighted to the central regions. Balmer lines have a strong 
centre-to-limb variation, becoming much weaker towards 
the limb, and therefore have a response to oscillations sim- 
ilar to that of radial velocities. Intensity measurements, 
on the other hand, give less spatial resolution because the 
centre-to-limb variation in intensity ('limb darkening') is 
relatively small. Bedding et al. went on to suggest that 
this difference in spatial response between simultaneous 
EW and intensity measurements could be used to deter- 
mine the £ value of a mode. 

They also investigated the temperature dependence 
of several other absorption lines (Fei, Can, G-band and 
Mg b feature) and found the Mg b feature and the Fe I lines 
to be especially useful for oscillation measurements. These 
spectral lines have much smaller centre-to-limb variations 
than the Balmer lines and so have spatial responses similar 
to that of intensity measurements. This means that mode 
identifications can be done with our observations by com- 
paring Balmer EW amplitudes with either intensity ampli- 
tudes or Fe I EW amplitudes. Modes with different degree 
I should occupy different locations in an amplitude-ratio 
diagram. 

A preliminary analysis of our observations was given 
by Viskum et al. (1997b), where only the ratio of Balmer- 
line EWs to intensity amplitudes was considered. In that 
paper, we found a clear grouping of the different modes in 
an amplitude-ratio diagram, and proposed a preliminary 
mode identification for the ten dominant frequencies. In 
this paper, we present a more thorough investigation of 
the different measurements, including the many Fei lines 
present in our spectra. 

Fig. H shows the amplitude ratios (Ha/Fel) versus 
(Ha/Intensity) for the eight dominant modes, obtained 
using the total Ha data set (Mt. Stromlo and La Silla). 
Note that amplitudes of oscillations in EW are measured 
in promille, while photometric amplitudes are measured in 
mmag. The data points come from Tables [| and |l while 
the 1-cr error bars are based on the tests on data subsets re- 
ferred to in Sec. || The intensity amplitudes are the mean 



of the amplitudes in the four Stromgren filters (uvby), in 
order to improve the signal-to-noise. Only the eight dom- 
inant modes are shown because of the large uncertainties 
on the remaining modes. 

The modes are clearly grouped into three regions. We 
identify the group at lower left as radial modes {£ = 0), fol- 
lowed by one dipole mode {1 = 1) and a group with £ > 2 . 
This trend follows the results of Bedding et al. ( 1996 ), 
who calculated that the amplitude ratio between different 
methods should become greater as £ increases. An impor- 
tant result from this diagram is that the Fe I lines indeed 
have spatial responses similar to that of intensity measure- 
ments, as predicted by Bedding et al. ( 1996 ), which shows 
that mode identification can be done from spectroscopic 
measurements alone. This is a very important result be- 
cause it allows one to obtain a mode identification based 
on a single data set, with a common window function and 
frequency resolution. 

One problem arises immediately: two of the modes 
in the lower-left group (f-2 — 280.42 ^Hz and J 3 = 
270.87 fiHz) are too closely spaced in frequency for both to 
be radial modes. Fig. ^| shows the same amplitude ratios, 
except that the Mt. Stromlo Ha data have been excluded. 
The separation of the modes is less clear. In particular, 
the ji mode has now moved away from the radial group 
towards the i = 1 group. We ascribe this to the change 
in resolution between the different data sets, which re- 
distributes the power between the two unresolved modes 
discussed in Sec. [| (/2 and /n). We have confirmed using 
simulations that, because of the presence of /n, the am- 
plitude ratio measured for the fi mode is very sensitive 
to the sampling and to their phase difference. We expect 
more reliable results when the Mt. Stromlo data are omit- 
ted because the spectroscopic data will then have a very 
similar window function to the photometric data. How- 
ever, our simulations show that even the remaining small 
differences in the sampling, which inevitably arise from 
the differences in statistical weights, will produce system- 
atic errors in the amplitude ratio. Of the modes identified 
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Fig. 4. Amplitude spectra in the four Stromgren filters 



by Breger et al. ( |1995| ), the pair f 2 and / n have by far the 
smallest separation (0.32 fiHz). Of those for which we have 
EW data, the mode involved in the next closest separation 
is /4, which is more than 2[iHz from /14. We therefore do 
not expect other modes in our figures to be affected to 
the same extent as f 2 , and we see that they are not. We 
conclude that the mostly likely identification is £ = for 
/ 3 and i = 1 for f 2 . 

Fig. [?] shows the amplitude ratios (H/3/Fel) versus 
(Ha/Fei). The diagram shows that the Ha and the H/3 
lines have essentially the same spatial response to the os- 
cillations. Together, the three amplitude-ratio diagrams 
indicate the sensitivity of different observables applied for 
mode identifications and the robustness of the method. 
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Fig. 5. Amplitude ratios (Ha/Intensity) versus (Ha/Fei) 
of the eight dominant modes in FG Vir. The EW ampli- 
tudes are calculated from the total data set (Mt. Stromlo 
and La Silla) . The intensity amplitudes are a mean of the 
four Stromgren filters (uvby). Labels show the mode fre- 
quencies in fiHz (truncated to integer values). 

Table 6. Mode identifications for FG Vir 

ID Frequency ^-ide ntifica tion 

fiYLz c/d Breger et al. (1995) This work 



h 
h 
h 
h 
h 
h 
h 
h 



147.18 
280.42 
270.87 
243.65 
229.95 
140.67 
111.76 
106.47 



12.72 
24.24 
23.41 
21.06 
19.87 
12.16 
9.66 
9.20 






2 
1 

2 

2 1 

2 1 



1 Possible g-mode. 

The mode identifications are summarized in Table |^ 
and Fig. pi In the table we also show the identifications 



given by Breger et al. (1995) and as seen the agreement is 
poor. This identification set is actually only one of several 
possible solutions selected by Breger et al. ( |1995| ) from a 
family of stellar models. They chose this solution, because 
it was in agreement with the constraints provided by the 
spectroscopic determination by Mantegazza et al. (1994), 
who found the fi mode to be radial. At the same time 



Breger et al. (1995) found from two-colour photometry a 
negative phase shift (<j>B-v — 4>v = —10.4 ± 2°) for the 
primary frequency (/1). This actually excludes the possi- 
bility of /1 being a radial mode, and indeed, this is also the 
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Fig. 7. Amplitude ratios (H/3/Fei) versus (Ha/Fei) of the 
eight dominant modes with la error bars. Only La Silla 
data 



result we obtain from the amplitude ratio diagrams. More- 



over, recent work by Breger et al. (1997) agrees with our 



identifications for all eight modes in Table |6j. The earlier 
disagreement clearly demonstrates the difficulty of classi- 
fying modes by using comparisons between stellar models 
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Fig. 8. Schematic diagram of the eight strongest modes in 
FG Vir, with our ^-identifications. The vertical axis shows 
the photometric amplitude in Stromgren y. 

and observed frequencies, without the help of additional 
information. 

Our identification of the strongest mode as dipolc 
(I = 1) is in good agreement with this mode's photo- 
metric negative phase difference obtained by Dawson et 
al. (1995). Another important result is that, as we show 
below, the fe mode is probably the fundamental radial 
mode (n = 1, £ = 0)). In this case, the two modes with 
lower frequencies (fa and /§) must be g-modes. These 
modes will be particular interesting for asteroseismology 
because they contain information about the convective 
core and convective overshooting. 

The above interpretation is based on the results ob- 



tained by Bedding et al. (1996), who considered the spe- 
cial case in which the rotation axis points towards the ob- 
server, so that contributions only arise from zonal modes 
(m = 0). Different orientations of the axis and the influ- 
ence of non-zero m values on the amplitude-ratio diagrams 
need to be investigated in detail. However, provided the 
dependence of amplitudes on m is smaller than the depen- 
dence on £, our identifications are likely to be correct. 

6. Comparison with a pulsation model 

6.1. Density of FG Vir 

Once the oscillation modes are identified, the next step 
in a seismic analysis is to compare the observed frequen- 
cies with models. Calculating a full theoretical model for 
FG Vir is beyond the scope of this paper. However, much 
progress can be made by scaling from existing model 
calculation s. He re we use calculations by Christensen- 
Dalsgaard ( 1993 ) for a 2.2M & star with solar metallicity, 
generated during investigations of the S Scuti variable k 2 
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Fig. 9. H-R diagram showing the evolutionary track of 
the 2.2 Mq model. The cross mar ks t he model most closely 
homologous to FG Vir (see Sec. 6.1). The vertical dotted 
lines show the 1-sigma range in T c g for FG Vir measured 
by Mantegazza et al (1994). The horizontal dotted lines 
show the 1-sigma range in luminosity determined from the 
Hipparcos distance (see Sec. |6~2| ). The pair of diagonal 
lines are contours of constant density and define the ±3er 
asteroseismic solution given in Table @. 



Bootis (Frandsen et al. 1995 ). Our aim is to make an ac 



curate estimate of the mean density of FG Vir using the 
radial modes. We do this by using a model that is homol- 
ogous to FG Vir. 

Two stars are homologous if their mass distributions 
are related by a simple spatial scaling factor (Kippen- 
hahn & Weigert 1990p , That is, suppose the two stars 
have masses M and M' and radii R and R' , and the mass 
enclosed by radius r (or r') is m (or m!). Then at all points 
where the relative radii are equal {r/R = r'/R'), the rela- 
tive enclosed masses will also be equal (m/M = m' /M 1 ). 
This being the case, and assuming the stars have the same 
chemical composition, the oscillation frequencies of one 
star can be found by multiplying the frequencies of the 
other star by a constant (equal to the square root of the 
ratio of their densities). 

The evolutionary track of the 2.2 M© model is shown 
in Fig. ^. The model is calculated from the zero-age main 
sequence (ZAMS) until the end of the core hydrogen- 
burning phase. Fig. |l0| shows the frequencies of the ra- 
dial modes of the 2.2M© model, plotted as a function of 
stellar density. As the model star evolves from right to 
left, the density and the frequencies both decrease. We 
only consider radial modes because those of higher degree 
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Fig. 10. Theoretical frequencies of the radial oscillation 
modes of a 2.2 M© star with solar metallicity (Z — 0.02), 
plotted as a function of stellar density (the star evolves 
from right to left). The dotted lines show the frequencies of 
the two observed modes in FG Vir that we have identified 
as radial. 
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Fig. 11. Frequency ratios for the radial modes of the 2.2 
M© model. The star evolves from right to left. The hori- 
zontal dotted line shows the ratio between the two modes 
in FG Vir that we have identified as radial. The vertical 
dotted lines shows the density of the corresponding model 
(0.203 po). That model should be homologous to FG Vir. 



have more complicated time dependencies (Christensen- 
Dalsgaard |1993| ; Frandsen et al. |1995| ). The horizontal dot- 
ted lines indicate the two observed modes in FG Vir that 
we have identified as radial (/3 and /e). Since the mass 
of FG Vir is somewhat less than 2.2 Mq, as we shall see, 
our aim is not to find an exact match to the model fre- 
quencies. Instead, we locate the stage in the evolution of 
the model at which the structure of the model is homol- 
ogous to that of FG Vir. This will occur when the model 
frequencies and those of FG Vir are related by a single 
multiplicative factor (which will be the square root of the 
density ratio) . The first step is therefore to identify a pair 
of radial modes in the model which match the observed 
frequency ratio {fa/ = 1.926). 
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Fig. 12. Scaled density for FG Vir versus model den- 
sity. For each radial mode, as a function of model density, 
we plot model density multiplied by (/obs //model) 2 - The 
model most closely homologous to FG Vir has a density 
of about 0.203 p© (vertical dotted line). The derived den- 
sity for FG Vir, shown by the horizontal dotted lines, is 
0.1645 ± 0.0005 p Q . 



To further establish the uncertainty in our density es- 
timate, we have constructed the plot in Fig. |l2] as fol- 
lows. For each of the three observed radial modes, we take 
each time step in the model evolution and use the ratio 
(/obs//modei) 2 to calculate the density that FG Vir would 
have if it were homologous to the model at that time step. 
This scaled density is plotted as a function of the model 
density. We see that the scaled density is fairly insensitive 
to exactly which stage in the model evolution is chosen. In 
other words, although the model most closely homologous 
to FG Vir has a density of 0.203 p© (vertical dotted line in 
Fig. |l2|), adjacent models do quite well (note the scales on 
the axes differ by more than a factor of 20). Indeed, the 
model with the same actual density as FG Vir has simi- 
lar absolute frequencies to FG Vir (to better than 0.5%), 
indicating that even this model is a reasonably good rep- 
resentation of FG Vir. 

On the basis of this discussion, we can set limits on 
the density of FG Vir as shown by the horizontal dotted 
lines in Fig. [l2[ which are chosen to span the three density 
estimates. We then have: 



Fig. [Tl| shows frequency ratios for the radial modes of P — 0.1645 ± 0.0005 p©. 



(3) 



the 2.2 M© model. We see that only one pair of modes 
matches the observed frequency ratio, which leads us to 
identify f§ with n = 1 and with n = A. The density of 
the corresponding 2.2 M© model is 0.203 p© and we select 
this as the model most closely homologous to FG Vir. To 
estimate the density of FG Vir, we simply multiply the 
density of the homologous model (0.203 p©) by the square 
of the ratio /obs//modei- The result is 0.164p©. 

We note in passing that the identification of f% = 
280.42 pHz as a radial mode (see Sec. ||) is inconsistent 



with the model frequencies. It is clear from Fig. 11 that 
the ratio /2//6 = 1-993 does not match the frequency ratio 
for any of the models. 

We can now attempt to identify other radial modes 
among the published frequencies. Based on the above den- 
sities, the radial mode with n — 2 s hould have a frequency 
of 185.8 pHz. Breger et al. ( 1998 ) observed a mode at 
186.04 pHz (/12) which is a good match. An exact fit to 
this observed frequency requires a density for FG Vir of 
0.1647p©. This differs from the above value by 0.3%, pro- 
viding an idea of the uncertainty in this process. It seems 
likely that /12 is a radial mode, but its amplitude is too 
weak for spectroscopic confirmation. 

Using the same argument, the radial mode with n = 3 
should have a frequency of 228.6 pHz. The two observed 
frequencies nearest this value are f§ — 229.95 p,Hz and 
f Q = 222.55 pHz (Breger et al. |1998| ). The best match is 
/5 but the discrepancy is 0.6% (1.2% in density) and our 
amplitude-ratio diagrams indicate that this mode is not 
radial. The frequency /g differs from the expected value 
by 2.7% and we have no information about its £ value. We 
therefore do not identify either of these as the radial n = 3 
mode. 



That is, assuming solar metallicity for FG Vir, we have 
determined the density to a precision of about 0.3%. It 
would clearly be desirable to obtain a spectroscopic mea- 
surement of the metallicity of this star, since this will affect 
the density at the level of about 1%. 

6.2. Distance and mass of FG Vir 

The effective temperature of FG Vir was de termined spec- 
troscopically by Mantegazza et al. ( 1994 ) to be 7500 ± 
100 K. From the model s use d in the seismic analysis of 
k 2 Boo (Frandsen et al. 1995| ), we have derived a relation- 
ship between luminosity, mass and effective temperature 
for stars with this effective temperature. For solar metal- 
licity (Z = 0.02) we find 



L 

Z©~ 



0.99 



M 

M© 



4.8 



(.. 



T, 



cil 



V 5777 K 



-0.8 



(4) 



Since L oc R 2 T^ S and mean density is proportional to 
M/R 3 , we then have: 



P ( Tpf! 

— ~ 1.0021 — — 
p Q ^5777 K 



R 

Rq 



-2.5833 



(5) 



Using the density given in Equation pi we can then es- 
timate the radius, luminosity and mass of FG Vir. Table 
shows the results for three values of effective temperature. 
The quoted uncertainties only include the uncertainty in 
the mean density. The values in Table |?] define a locus of 
constant density in the H-R diagram, as shown in Fig. |[ 
The portion of this curve that lies between the extremes 
of valid temperatures defines the solution. 
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Table 7. Parameters of FG Vir for three different values 
of temperature 





7400 K 


7500 K 


7600 K 


R/Rq 


2.215 ± 0.003 


2.227 ± 0.003 


2.238 ± 0.003 


L/Lq 


13.21 ± 0.03 


14.08 ± 0.03 


15.00 ± 0.03 


M/Mq 


1.788 ± 0.008 


1.816 ± 0.008 


1.844 ± 0.008 


logs" 


3.999 ± 0.001 


4.002 ± 0.001 


4.004 ± 0.001 



a units of g are cms 



Table 8. Parameters of FG Vir based on asteroseismic 



density 12 


p/pe 


0.1645 ± 0.0005 


R/Rq 


2.227 ± 0.012 


l/l q 


14.1 ± 0.9 


m/m q 


1.82 ± 0.03 


logs 


4.002 ± 0.003 


distance (pc) b 


84 ± 3 



a assuming T cB = 7500 ± 100 K and Z = 0.02. 

b using L/Lq, V = 6.57, BC = -0.07 and M bo i, = 4.75. 



Table 9. 

distance" 



Parameters of FG Vir based on Hipparcos 



distance (pc) 
L/L Q b 

R/Rq 

M/Mq 

p/pQ 

log g 



83 ± 5 
13.8 ± 1.6 
2.20 ± 0.15 
1.82 ± 0.03 
0.17 ± 0.03 
4.0 ± 0.05 



a assuming T cff = 7500 ± 100 K and Z = 0.02. 
b using distance, V = 6.57, BC = -0.07 and M\ 



bol,0 



: 4.75. 



The main sources of uncertainty in this process arise 
from imprecise knowledge of the temperature and metal- 
licity and from the use of the relation between L, M and 
T c g (Equations |] and ||). Assuming that the main uncer- 
tainty is the temperature, we find the properties of FG Vir 
listed in Table ||. The errors are not independent: increas- 
ing the temperature will increase the luminosity, mass and 
radius. 

The derived distance is in excellent agreement with the 
Hipparcos value of 83 ± 5 pc, as is also clear from Fig. [| It 
is important to note that if the temperature had been bet- 
ter known, e.g. to ±30 K, we would have obtained an error 
on the distance of only ±1 pc (internal error). The surface 
gravity agrees with the value of logg = 3.9 ± 0^ deter- 
mined spectroscopically by Mantegazza et al. ( 1994 ). The 
high precision of our log g determination arises because g 
is given by pR. 

Table ^| shows the fundamental properties for FG Vir 
obtained by combining the Hipparcos distance with Equa- 
tions H and H, without using the asteroseismically deter- 
mined density. It can be seen that the asteroseismic prop- 



erties for FG Vir are far more precise than those obtained 
by using the Hipparcos data, even with 100 K uncertainty 
in the temperature. Such calculations show clearly the po- 
tential of asteroseismology. The present analysis is based 
on three £ = modes out of more than 20 known fre- 
quencies. We have identified modes with I — 1 and £ = 2 
and these modes will contain accurate information on the 
detailed evolutionary stage of FG Vir. This should pro- 
vide strong constraints on models computed specifically 
for FG Vir, resulting in a fine-tuned model. The effect 
from metallicity should also be investigated. 

To check our identification of I = modes, we com- 



pared all 24 frequencies published by Breger et al. fll998| ) 
with the frequency ratios of radial modes in the 2.2M Q 
model. The highest correlation occurred with the follow- 
ing solution: f 6 (n = 1), f 12 (n = 2), f 5 (n = 3) and / 3 
(n = 4), with a scatter between observed and best-fit the- 
oretical frequencies of about 1 /iHz. These identifications 
are in full agreement with our own results and give us con- 
fidence that we have correctly identified the radial modes 
in FG Vir. 

Finally, we note that the strongest mode (/i) is likely 
to be n = 1, 1 = 1 and is probably undergoing an avoided 
crossing with a g-mode. This would make this mode very 
interesting for testing a more detailed model. 



7. Rotation period of FG Vir 

A periodic trend was evident in the raw EW data for all 
metal lines and in the two Balmcr lines. The same peri- 
odicity was seen in all four photometric light curves. The 
variation is clearly seen after removal of the known pulsa- 
tion signal and smoothing (Fig. |l3|) . 

There appears to be a variability with a period of 
3.5 ± 0.2 days and a decaying amplitude. Note that the 
Ha time series is composed of data from telescopes at two 
different sites, while the Fei and Ha time series were ob- 
tained using different reduction programs. Furthermore, 
the photometric time series were obtained with a third 
telescope. All time series are in phase, indicating that the 
variation is caused by the same phenomenon. 

We propose that this 3.5-day variation is intrinsic to 
the star and reflects the rotation period of FG Vir. Using 
the asteroseismic radius in Table ||, we then obtain v eq = 
33±2 km s _1 for the equatorial velocity. Using the value of 
vsmi = 21 ± lkms -1 found by Mantegazza et al. (1994) 
gives i = 40° ± 4° for the inclination angle, in agreement 
with the value of i — 31° ± 5° found by Mantegazza et al. 
from a moments analysis of /s and fr. 



8. Summary of results 

We have investigated a new technique to measure the os- 
cillations in S Scuti stars via changes in the equivalent 
widths of absorption lines. An important advantage of 
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Fig. 13. The raw time series for (a) the photometric light 
curve in the u filter; (b) the Ha line EW (small symbols: 
Mt. Stromlo data, large symbols: La Silla data); (c) the 
combined Fel EW time series. The dotted curve shows a 
damped oscillation with a period of 3.5 days and with a 
mean lifetime (the inverse of the damping constant) of 14 
days. 



this new technique is that only medium-dispersion spec- 
tra are needed, which makes the method suitable for small 
and medium-sized telescopes and for multi-site campaigns. 
Our main results are summarized below. 

— Our detection of oscillations in FG Vir from 
equivalent-width measurements of Ha, H/3 and Fel 
lines is an important confirmation of the method devel- 



oped by Kjeldsen et al. (1995) to search for solar-like 
oscillations. 

— From the ratios between oscillation amplitudes mea- 
sured in EW and the four Stromgren filters (uvby), we 
have identified ^-values for eight modes in FG Vir. 

— We suggest that the two lowest-frequency modes are g- 
modes, while the strongest mode (147.2/iHz) is a dipole 
mode. 

— By comparing the frequencies of radial modes with 
model calculations, we obtained a precise density and 



derived a distance that is in excellent agreement with 
the Hipparcos value. 
— We detected a long-period variation in the time series 
with a period of 3.5 days, which we propose is caused 
by rotation of the star. 

Acknowledgements. We are grateful to M. Breger for provid- 
ing us with oscillation data of FG Vir prior to publication 
and for useful discussions. We are grateful to J.V. Clausen for 
his assistance during the observations with the Danish 0.5m 
telescope. We thank F. Kupka and M. Gelbmann for provid- 
ing us with spectral line identifications and synthetic spectra 
for FG Vir, and J.G. Robertson for many useful comments on 
the manuscript. This work was supported in part by the Dan- 
ish National Research Foundation through its establishment 
of the Theoretical Astrophysics Center, and also by the Aus- 
tralian Research Council. We are grateful to the HIPPARCOS 
group for making their catalogue available. This research has 
made use of the Simbad database, operated at CDS, Stras- 
bourg, France. 



References 

Bedding, T.R., Kjeldsen, H. 1998, in: Proceedings of the Tenth 
Cambridge Workshop on Cool Stars, Stellar Systems and 
the Sun, R.A. Donahue & J. A. Bookbinder (eds.), Cam- 
bridge, USA, in press 

Bedding, T.R., Kjeldsen, H., Reetz, J., Barbuy, B. 1996, MN- 
RAS, 280, 1155 

Breger, M., 1997, private communication and Introduction to 
Delta Scuti Star Newsletter, M. Breger (ed.), Vienna 1997, 
Vol. 11, p. 3 



( tittp : / /dsn. ast .univie . ac . at/dssn/nlll/chapter_L .html ) 



Breger, M., Handler, G., Nather, R.E., Winget, D.E., et al. 

1995, A&A, 297, 473 
Breger, M., Zima, W., Handler, G., Poretti, E., et al. 1998, 

A&A, 331, 271 

Christensen-Dalsgaard, J. 1993, in: Inside The Stars, W.W 
Weiss & A. Baglin (eds.), ASP Conf. Ser. 40, 483 

Dawson, D.W., Breger, M., Lopez de Coca, P. 1995, PASP, 
107, 517 

Frandsen, S., Jones, A., Kjeldsen, H, Viskum, M., Hjorth, J., 
Andersen, N.H., Thomsen, B. 1995, A&A, 301, 123 

Frandsen, S., Balona, L.A., Viskum, M., Koen, C, Kjeldsen, 
H. 1996, A&A, 308, 132 

Gelbmann, M. 1995, Comm. Asteroseismology 89 

Gelbmann, M., Kupka, F., Weiss, W.W., & Mathys, G. 1997, 
A&A, 319, 630 

Gray, D.F. 1992, Cambridge Astrophys. Ser. Vol. 20, The ob- 
servation and analysis of stellar photospheres, 2nd edn. 
Cambridge Univ. Press, Cambridge 

Guzik, J. A. & Bradley, P. A. 1995, Baltic Astronomy, 4, 442 

Kippenhahn, R. & Weigert, A. 1990, Stellar Structure and Evo- 
lution, Springer- Verlag, Berlin 

Kjeldsen, H, Bedding, T.R., Viskum, M., Frandsen, S. 1995, 
AJ, 109, 1313 

Mantegazza, L., Poretti, E., Bossi, M. 1994, A&A, 287, 95 
Press, W.H., Flannery, B.P., Teukolsky, S.A., Vetterling, W.T. 
1989, Numerical Recipes, Cambridge Univ. Press, Cam- 
bridge 



M. Viskum et al.: The 8 Scuti star FG Vir 



13 



Viskum, M., Baldry, I.K., Kjeldsen, H., Frandsen, S., Bed- 
ding, T.R. 1997a, in: Sounding Solar and Stellar Interiors, 
J. Provost & F.-X. Schmieder (eds), IAU Symp. no. 181, 
Kluwer, Dordrecht, Holland, in press 

Viskum, M., Dall T.H., Kjeldsen, H., Baldry I.K., Bruntt, H., 
Frandsen, S., Bedding, T.R. 1997b, in: A half century of 
stellar pulsation interpretations: A tribute to Arthur N. 
Cox, J. A. Guzik & P. Bradley (eds), ASP Conf. Ser. 135, 
465 



